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ABSTRACT 


A unified approach for solving incompressible flows has been investigated in this study. 
The numerical CTVD (Centered Total Variation Diminishing) scheme used in this study 
was successfully developed by Sanders and Li for compressible flows, especially for the 
high speed. The CTVD scheme possesses better mathematical properties to damp out 
the spurious oscillations while providing high-order accuracy for high speed flows. It 
leads us to believe that the CTVD scheme can equally well apply to solve incompressible 
flows. Because of the mathematical difference between the governing equations for 
incompressible and compressible flows, the scheme can not directly apply to the 
incompressible flows. However, if one can modify the continuity equation for 
incompressible flows by introducing pseudo-compressibility, the governing equations for 
incompressible flows would have the same mathematical characters as compressible 
flows. The application of the algorithm. to incompressible flows thus becomes feasible. 

In this study, the governing equations for incompressible flows comprise continuity 
equation and momentum equations. The continuity equation is modified by adding a 
time-derivative of the pressure term containing the artificial compressibility. The modified 
continuity equation together with the unsteady momentum equations forms a hyperbolic- 
parabolic type of time-dependent system of equations. Thus, the CTVD schemes can be 
implemented. In addition, the physical and numerical boundary conditions are properly 
implemented by the characteristic boundary conditions. 

Accordingly, a CFD code has been developed for this research and is currently under 
testing. Flow past a circular cylinder was chosen for numerical experiments to determine 
the accuracy and efficiency of the code. The code has shown some promising results. 
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INTRODUCTION 


GOVERNING EQUATIONS 

The two-dimensional incompressible-flow equations expressed in conservative variables, 
p, pu, and pv are as follows: 


dx dy 


(D 


dptj + 9pu 2 + 8p w = _dp , p, , (fyu + <Ppu_ \ (2) 

dt dx dy dx p dx 2 dy 2 


dpv + dpuv + Bpv 2 = _ 9p , p, , &pv + 8 2 pi^ (3) 

dt dx dy dy P dx 2 dy 2 


where u and v are velocity components in the x and y directions, respectively. P is the 
static pressure, p denotes the dynamic viscosity, p is the fluid density. The 
incompressible governing equations (1) through (3) are mathematically classified as 
elliptic partial differential equations while compressible governing equations are hyperbolic 
partial differential equations. Because of the mathematical difference between the 
hyperbolic and elliptic partial differential equations, a well-developed numerical schemes 
for compressible flows can not apply directly to solve incompressible flows. However, if 
one modifies the continuity equation given in equation (1) by introducing artificial 
compressibility, the resulting incompressible governing equations are of hyperbolic type. 
The modified continuity equation including a time-derivative of the pressure term 
containing the artificial compressibility is given in equation (4). 


3p + appiy + 9ppv =0 (4) 

dt dx dy 
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where p is known as the pseudo-compressibility constant and has the dimension as 
velocity square. 


Equations (2) through (4) are transformed onto generalized curvilinear coordinates, \ and 
ri given by 




(5) 


r\=r\{x,y,t) 


(6) 


The governing equations are then given by 

9Q x a(E-E (7) 

dx dr) 


where 



P 

ppu 

P9 V . 


(8) 


E= 


P uU+\j) 

[pW^ 


(9) 



Mv-n,) 

pUV+T\,p 
P VV+T[yP_ 


( 10 ) 
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( 11 ) 


0 






o 

(^ + ^y)^ + (nx + n5)^„ 

tox + W^+(Tl!|+Tl5)K, 


( 12 ) 


U and V are contravariant velocities. J is the Jacobian of the transformation, the metrics 
of the transformation are 




H f s dl 
dx y ~ dy 


and 


_5tl 

dx' 


_3ti 




(13) 


(14) 


U and V are contravariant velocities and given as follows: 

U-ljJ+l y v, V=r\^+r\ y v (1 5) 


The transformed governing equation given in equation (7) is then solved by a finite 
volume method. Accordingly, the conservation laws given below in integral form is used 
in the formulation. 


— fQdV+ j' DdA=0 (16) 
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D is derived from equation (7) and a second-order flux tensor defined in % and rj. 


NUMERICAL SCHEME 

A CTVD scheme developed by Sanders and Li (Ref. 1) is adopted to determine the 
interfacial variables for the finite volume formulation. The CTVD differencing has 
distinctively desired properties in overcoming the spurious oscillations and odd- and even- 
point decoupling in the solution which are caused by the use of central differencing. This 
scheme has demonstrated the superiority in reducing CPU time while providing up to 
fourth order accuracy. In addition, the implementation of these algorithms is simple and 
no tuning parameters are needed. 

By applying the interfacial variables derived from the CTVD scheme to the conservation 
laws given in equation (16) numerical fluxes are evaluated on the cell faces based on the 
spatial discretization. Equation (16) is further reduced to a system of time dependent 
differential equations by the method of lines. The system of differential equations can be 
integrated by a number of relaxation methods. In this study, the explicit Runge-Kutta 
method is employed for solutions. 

The implementation of wall and far-field boundary conditions is crucial to the numerical 
scheme. The flows are treated as subsonic. Therefore, the proper far-field boundary 
conditions must be imposed. Characteristic boundary conditions are used and integrated 
into CTVD scheme to evaluate fluxes. 

NUMERICAL EXPERIMENTS 

A CFD code based on the CTVD scheme has been developed to solve two-dimensional 
incompressible flows. Numerical testings have been conducted to solve flow past a 
circular cylinder problem. The far-field boundary conditions are imposed at a distance 1 0 
times of radius from the center of the cylinder. The initial conditions are given as follows: 

p =2.0 Ib/ft 2 
U 0 =1 .0 ft/sec. 
p=1 .0 slug/ft 3 

R=1.0 ft., is the radius of the cylinder. J3=5, 10,30,40, and 50 ft 2 /sec 2 are utilized in 
numerical testings. The O-grid system is used to discritize the flow domain. All numerical 
testings are based on a grid of 61x82 nodes. A periodic boundary condition is imposed 
around the branch cut. The maximum CFL number can be used in each testing is 0.6. 



According to the numerical results, [3 value plays very important role in the numerical 
calculation. When p is greater than 1 0 and less than 40, the solution converges. A typical 
second-order pressure field is given in Figure 1 after 500 iterations. 

The numerical testings have not finished yet. More testings are need to investigate 
properly imposing wall boundary conditions. How to increase CFL number in the 
calculation is also a subject for more research. 

CONCLUDING REMARKS 

A research CFD code has been developed for incompressible flows using the CTVD 
scheme. The numerical results have shown that the CTVD scheme can apply to 
incompressible flows by the use of a modified continuity equation. A converged solution 
can be obtained by choosing a proper pseudo-compressibility constant. A further 
investigation on the CFL number is necessary in order to increase the efficiency of the 
method. The implementation of diagonal ADI method to solve the reduced set of ordinary 
differential equations would increase the efficiency of the method. 
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